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We introduce a set of four new publicly available A''-body simulations, the most recent addi- 
tions to the Texas P^M Database. Our models probe the less studied parameter space region of 
moderate volume (100/i~^Mpc box) combined with fine mass resolution (oc 10^^ Mq, roughly 
comparable to a galaxy), making these simulations especially suitable for study of major 
large-scale structure (LSS) features such as voids, and for comparison with the largest three- 
dimensional redshift surveys currently available. Our cosmological models (LCDM, TOCDM, 
OCDM, TCDM) are all C05£-normalized, and when possible (LCDM and TOCDM) also 
cluster-normalized, based on the X-ray cluster M-T relation. The COBE- and cluster- 
normalized LCDM model reiterates the attractiveness of this currently favored model which 
does not require the introduction of tilt in order to fit the constraints imposed by observations 
of other cosmological parameters. 



1 Introduction 

A^-body simulations are an essential tool for probing LSS and galaxy formation. As large, high- 
resolution simulations are computationally costly, one has to carefully consider the added effort 
resulting from increasing the simulations' volume or from improving their resolution. In this 
context, most simulations gravitate towards a design stressing either of these two conflicting 
goals. The largest three-dimensional redshift surveys currently available are situated somewhere 
in between these two extremes: a M^, galaxy in the CfA2 survey is visible out to 100/i~^Mpc. A 
simulation designed to match these surveys must have both the required resolution to identify 
the dark matter (DM) halos associated with such galaxies, and this moderately large volume. 

The other essential consideration of simulation design is the choice of cosmological models 
probed. Ideally, one would want to examine a certain range of the relevant cosmological param- 
eters {Hq, ^0) ^^BO) "^s)) but this is often not an attainable goal. In this work we focus 
on cosmological models with currently favored values of Hq, fig (and Aq), and require that all 
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models will be COi?£^-normalized. The above constraints can still be fulfilled through a variety 
of primordial power spectrum tilt n and erg combinations. We attempt to achieve also cluster 
normalization, hence determining the value of (and thus fixing a tilt value). 

2 Models 

These simulations were designed with the goal of maximizing their volume while still being able 
to resolve DM halos associated with galaxies. Adopting a mass-to-light ratio of IOOMq/Lq, 
we thus require that we will be able to identify Afhaio ^ IO^^Mq halos. A simulation with 140^ 
particles in a 100/i~^Mpc box will have Mparticie = 1-01 x IO^^JIq /^^^M©- Identifying all DM 
halos with 20 particles or more, this parameter specification matches the stated requirements. 

Failure to identify halos all the way down to this resolution renders such simulations inap- 
propriate for the study of many LSS features. Specifically, it would be impossible to identify 
correctly voids in such simulations. Such studies often originate from DM simulations, and use 
some form of halo populating scheme in order to match the observed properties of the distribution 
of galaxies. These populating schemes (also known as bias recipes) assign a number of galaxies 
to each halo. But for such schemes to have a fighting chance at reproducing the distribution of 
galaxies, one must initially know the locations of the halos that should be populated — including 
halos that would be populated by just one galaxy. 

We estimated that 20 particles per halo is the minimal number required in order to reliably 
identify halos. As this is getting close to the simulation's mass resolution limit, we tested the 
lower end of our DM halo mass function by constructing a matching set of smaller simulation 
boxes with Mparticie Order of magnitude smaller □ (see Table I). We can then compare the 
number density of DM halos in the two sets of boxes and see if indeed we manage to recover the 
correct number of small halos in the larger simulation boxes. 

All simulations were started at an initial redshift Zi = 24 and evolved over 600-1000 timesteps 
using a P^M code over a 256^ grid. More details on the simulations can be found elsewhere. U'H 

In Table || we summarize the cosmological parameters defining our models. Column 1 
indicates the models' acronyms. Columns 2-5 detail the models' values of the Hubble parameter 
/i, density parameter cosmological constant Aq, and tilt n. Column 6 lists the theoretical 
^cluster ]-,g^gg(^ X-ray cluster temperatures. Column 7 lists C7g°°*, the actual value corresponding 
to each cosmological model. If these two values match, we state that the model is cluster- 
normalized (Column 8). 

We fixed all models to be COiJ-B-normalized with Tcmb = 2.7i^. assuming no contribution 
from tensor modes. Also, we used Obo = 0.015/i~^ throughout, □ in concordance with th£ 
Texas P^M Database. When practical (all models but TCDM) we adopt a Hubble constants 
Hq = 65kms~^Mpc~^. In addition to COBE normalization, we have also tried to achieve 
cluster normalization. Using the X-ray cluster M-T relation Q we derived the required erg value 
for each of our models. For each combination of Oq, Aq, and ffo; we computed the tilt required 
in order to achieve cluster normalization and examined whether it is acceptable in view of the 
limits allowed by the 4-year COBE data. 

For two of the models — LCDM and TOCDM — we found an acceptable tilt value and man- 



Table 2: Model Parameters 
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Figure 1: CMB Angular Power Spectra: Observations vs. Models 



aged to achieve cluster normalization. The two other models are not cluster-normalized. The 
OCDM model was designed as a direct companion to the LCDM model, where all the cosmolog- 
ical parameters in both these models — except the value of Aq — are the same. The TCDM model 
is our best attempt with an = 1 model, where we used the lowest possible fig value which 
does not require h < 0.55 or n < 0.7. For comparison, we have also included in Table |2| (and in 
Fig. H) the familiar SCDM model, although it is not one of the cosmological models simulated. 

3 Results 

In Fig. |l] we compare the theoretical CMB angular power spectra of the models simulated here 
with the recent BOOMERanCi anisotropy measurements. In Fig. |2| we present cumulative halo 
mass functions for the two sets of models simulated. Our two cluster-normalized models. LCDM 
and TOCDM, reproduce similar mass functions. The observational point in the figurelll n{T > 
4.0keV) = 1.5±0.4xl0~*^/i3Mpc~2 is in good agreement with the LCDM cluster abundance. The 
TOCDM curve follows closely the LCDM curve, but for the former cosmology the observational 
point would be shifted along the horizontal axis by a factor 0.3^/^. However, it should be noted 
that there are still significant uncertainties associated with both the observational measurements 
of cluster abundance and the theoretical modeling of the M~T relation, i 

There are two curves for each cosmological model in Fig. |2| — one representing the mass 
function as measured in the 100/i~^Mpc box, the other measured in the 40/i~^Mpc box. As 
illustrated in the figure, for each model there is excellent agreement between the two curves. 




Figure 2: Halo Mass Functions 



4 Summary 

In this paper we introduce two matching sets of four cosmological models. We derive halo mass 
functions for all models and use the small box, high resolution simulations in order to verify 
the validity of the mass function in the large box for halos as small as ~ W^'^Qq h~^MQ. The 
simulations presented here are unique as they both cover a volume comparable to current large 
three-dimensional redshift surveys and at the same time resolve cluster masses down to M^. 
While the simulations were designed mostly in order to achieve COBE — and, when possible 
(LCDM and TOCDM), also cluster — normalization, they also serve to demonstrate the attrac- 
tiveness of the LCDM model. Requests for the simulations presented in this paper, or for other 
data from the Texas P^M Database, should be sent to database@galileo.as.utexas.edu. 
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